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Abstract: Interactions of virtual photons are analyzed in terms of photon struc- 
ture. It is argued that the concept of parton distribution functions is phenomeno- 
logically very useful even for highly virtual photons involved in hard collisions. The 
role of longitudinal photons for proper interpretation of jet cross-sections in the 
region of moderate virtualities accessible at HERA is investigated. 
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'. 1 Introduction 

o : 

In quantum field theory it is difficult to distinguish effects of the "structure" from those of 
"interactions". Within the Standard Model it makes good sense to distinguish fundamental 
Q^; particles, which correspond to fields in its lagrangian £ SM (leptons, quarks, gauge and Higgs 
bosons) from composite particles, like atoms or hadrons, which appear in the mass spectrum 
but have no corresponding field in £sm- For the latter the use of parton distribution functions 
(PDF) to describe their "structure" appears natural, but the concept of PDF turns out to be 
.£h ! phenomenologically useful also for some fundamental particles, in particular the photon. PDF 
are indispensable for the real photon due to strong interactions between the pair of quarks 
to which it couples electromagnetically. For massless quarks this coupling leads to parallel 
singularities, which must be absorbed in PDF of the photon. For nonzero photon virtualities 
there is no true singularity associated with the coupling 7 — > qq and therefore no real need for 
PDF. The main aim of this paper is to advocate the use of PDF even for the virtual photon. 

2 PDF of the real photon 

The factorization scale dependence of PDF of the real photon is determined by the system of 
coupled inhomogeneous evolution equations 

d£(M 2 ) 



dlnM 2 
dG(M 2 ) 
din M 2 
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k q + P qq ® £ + P qG ® G, (1) 
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1 Contribution to the Workshop Monte Carlo generators for HERA Physics, DESY, 1998-1999 
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Figure 1: Diagrams defining the pointlike parts of quark and gluon distribution functions in LL 
approximation. The resummation involves integration over parton virtualities Mq <t< M 2 . 



for singlet, nonsinglet and gluon distribution functions. The spliting functions k q , k G and Py 
admit expansions in powers of a s , 
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general solution of these equations can be written as the sum of a particular solution of the 
full inhomogeneous equations and the general solution of the corresponding homogeneous ones, 
called hadronic (or VDM) part. A subset of solutions of the inhomogeneous evolution equations 
resulting from the resummation of diagrams in Fig. [3] defines the so called pointlike (PL) 
parts. This resummation softens the x— dependence of q PL (x, M 2 ) with respect to the term 
(a/2ir)k^l(x), corresponding to the simple 7 — > qq splitting. A general solution of eqs. ([1|-|2|) 
can thus be written as (D = q, q, G) 



D(x,M 2 



D PL (x,M 2 ) + D 



VDM, 



x,M 2 



D = q,q,G. 



(7) 



The important point to keep in mind is the fact that there is an infinite number of pointlike 
solutions q PL (x, M 2 ), G PL (x, M 2 ), differing merely by the initial scale M at which they vanish. 
Consequently, the separation of quark and gluon distribution functions into their pointlike and 
VDM parts is ambiguous and these concepts have thus separately no physical meaning 0. 



2.1 Properties of Schuler— Sjostrand parametrizations 

Practical aspects of the ambiguity in separating PDF into their VDM and pointlike parts can 
be illustrated on the properties of SaSID and SaS2D parametrizations 0,0, which differ by 
the choice of initial M : M = 0.6 GeV for SaSID, M = 2 GeV for SaS2D. What makes 
the SaS approach particularly useful for our discussion is the fact that it provides separate 
parametrizations of the VDM and pointlike parts of both quark and gluon distributions. In 

2 For brevity the terms "pointlike quarks" and "pointlike gluons" will hence be employed to denote pointlikc 
parts of quark and gluon distribution functions of the photon. 
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Figure 2: The u and c quark and gluon distribution functions of the real photon for SaSID 
(upper solid curves) and SaS2D (upper dotted curves) parametrizations at M 2 = 100 GeV 2 . The 
VDM and pointlike ( solid and dotted curves peaking at large x ) parts of both parametrizations 
are plotted separately. For quarks the splitting terms corresponding to SaSID (upper dashed 
curves) and SaS2D (lower dashed curves) parametrizations are overlayed. 



Fig. we compare distribution functions xu(x,M 2 ), xc(x,M 2 ) and xG(x, M 2 ) as given by 
SaSID and SaS2D parametrizations at M 2 = 100 GeV 2 . To see the effects of the resummation 
of multiple parton emission we also plot the corresponding splitting terms 

<f&{x, Ml M 2 ) ee f 3e 2 (x 2 + (1 - x) 2 ) In ^. (8) 

Z7T M 

In Fig. § we illustrate in two ways the scale dependence of VDM and pointlike parts of the 
same three distribution functions. In the upper six plots we compare them as a function of x 
at M 2 = 25, 100, 1000 GeV 2 , while in the lower three plots the same distributions are rescaled 
by the factor ln(M 2 /M 2 ). Finally, in Figs. | and | we compare LO SaS predictions for two 
physical quantities: F% and effective parton distribution function D e g, relevant for approximate 
calculations of jet production 

n f 

Fj(x,Q 2 ) = ^2xe 2 % (*,Q 2 ), (9) 
i=i 

D cS (x,M 2 ) = ^(g,(x,M 2 ) + g l (x,M 2 )) + ^G(x,M 2 ). (10) 
i=i 

Figures |2|-|5] illustrate several important properties of PDF of the real photon: 



• There is a huge difference between the relative importance of VDM components in SaSID 
and SaS2D parametrizations: for SaS2D the VDM components of xu(x, M) and xG(x, M) 
are dominant up to x = 0.75, whereas for SaSID the pointlike one dominates already 
above x = 0.1! 

• Factorization scale dependence of VDM and pointlike parts differ substantially. VDM 
components exhibit the pattern of scaling violations typical for hadrons, whereas the 
pointlike parts of both quark and gluon distribution functions rise with M for all x. For 
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Figure 3: Factorization scale dependence of parton distributions functions u(x, M), c(x, M) and 
G(x, M) of the real photon. Dashed curves correspond, in the order indicated by the arrows, to 
pointlike and solid to V DM parts of these distributions at M 2 = 25, 100 and 1000 GeV 2 . The 
meaning of arrows is the same for all parton distribution functions. In the lower part SaSID 
quark distribution functions xu(x, M 2 , P 2 ) scaled by 1ii(M 2 /Mq) for the real photon (left) and 
by \n(M 2 /P 2 ) for the virtual one, are plotted and compared to the predictions of the splitting 
term jj^). 
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Figure 4: F^x^Q 2 ) as a function of x for Q 2 = 25, 100, 1000 GeV 2 as given by SaSID (solid 
curves) and SaSD2 (dotted curves) parametrizations. The full results are given by the upper, 
the VDM and pointlike contributions parts by two lower curves. The dashed curves describe 
the contributions of the splitting term ^) with M = 0.6 GeV for SaSID and M = 2 GeV for 
SaS2D. 
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Figure 5: D e fj(x,M 2 ) as a function of x for M 2 = 25, 100, 1000 GeV 2 for the real photon as 
given by SaSID and SaSD2 parametrizations. Solid curves show the full results, dashed ones 
the VDM and dotted ones the pointlike parts. The pointlike part is further separated into the 
contributions of pointlike quarks and pointlike gluons, denoted by two dash-dotted curves (those 
peaking at low x correspond to gluons). Thick dashed curves correspond to the splitting term 
(^). The thin solid curves in SaS2D plots show the full results of the SaSID parametrization. 
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pointlike gluons this holds despite the fact that G PL (x, M 2 ) satisfies at the LO standard 
homogeneous evolution equation and is due to the fact that its evolution of is driven by 
the rise of E PL (x,M 2 ). 

• As the factorization scale M increases the VDM parts of both quark and gluon distribution 
functions decrease relative to the pointlike ones, except in the region of very small x. 

• Despite huge differences between SaSID and SaS2D parametrizations in the decomposi- 
tion of quark and gluon distributions into their VDM and pointlike parts, their predictions 
for physical quantities F% and D efi are quite close. 

• The most prominent effect of multiple parton emission on physical quantities appears to 
be the contribution of pointlike gluons to jet cross-sections in the region x 7 < 0.5. 

3 PDF of the virtual photon: do we really need them? 

For sufficiently virtual photon the initial state singularity resulting from the splitting 7* — > qq 
is shielded off by the nonzero initial photon virtuality P 2 and therefore in principle the concept 
of PDF does not have to be introduced. Nevertheless, even in such circumstances PDF turn 
out to be phenomenologically very useful because their pointlike parts include the resummation 
of parts of higher order QCD corrections, and the VDM parts, though decreasing rapidly with 
increasing P 2 , are still dominant at very small x 1 . Both of these aspects define the "nontrivial" 
structure of the virtual photon in the sense that they are not included in the splitting term (|S|) 
and thus are not part of existing NLO unsubtracted direct photon calculations. 

3.1 Virtuality dependent PDF 

In QCD the nonperturbative effects connected with the confinement are expected to determine 
the long-range structure of the photon and hence also the transition between the virtual and 
real photon. As for the real photon, we recall basic features of SaS parametrizations of PDF of 
the virtual photon. We refer the reader to the web version of this contribution for the figures 
illustrating the following observations. 

• Both VDM and pointlike parts drop with increasing P 2 , but VDM parts drop much faster. 

• With increasing P 2 the realative importance of VDM parts of both quark and gluon 
distribution functions decreases rapidly. For M 2 > 25 GeV 2 the VDM parts of both 
SaSID and SaS2D parametrizations become practically negligible already at P 2 w 3 
GeV 2 , except in the region of very small x < 0.01. Hence, also the ambiguity in the 
separation (0) is largely irrelevant in this region. 

• The general pattern of scaling violations remains the same as for the real photon, except 
for a subtle difference, best visible (see Fig. [3]) when comparing the rescaled PDF for 
P 2 = with those at P 2 = 1,3 GeV 2 . While for P 2 = the SaSDl parametrizations of 
quark distribution functions soften with increasing M 2 and intersect the splitting term at 
x ~ 0.65, for P 2 > 1 GeV 2 they rise with M 2 for all x and stay above the splitting term. 
These properties reflect the fact that SaS parametrizations of PDF of the virtual photon 
do not satisfy the same evolution equations as PDF of the real one. 
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• Pointlike quarks dominate D e g(x, P 2 , M 2 ) at large x, while for x < 0.5, most of the 
pointlike contribution comes from the pointlike gluons. In particular, the excess of the 
pointlike contributions to D e g over the contribution of the splitting term, observed at 
x < 0.5, comes almost entirely from the pointlike gluons! 

• For x > 0.6 the full results are below those given by the splitting term (|8|) with Mq = P 2 
and one therefore expects the sum of subtracted direct and resolved contributions to jet 
cross-sections to be smaller than the results of unsubtracted direct calculations. 

Jet production in ep collisions at HERA in the region of photon virtualities P 2 > 1 thus offers a 
promising opportunity for the identification of nontrivial aspects of PDF of virtual photons at 
both small (but not very small) and large values of x. All these conclusions were obtained within 
the simple LO formalism, but as shown in the next Section they hold within the framework of 
NLO parton level calculations as well. 

3.2 Should we care about 7£? 

Most of the existing phenomenological analyses of the properties and interactions of virtual 
photons as well as all available parametrizations of their PDF concern transverse photons 
only. Neglecting longitudinal photons is a good approximation for y — > 1, where the flux 
f2{y,P 2 ) — > 0, as well as for very small virtualities P 2 , where PDF of 72 vanish by gauge 
invariance. But how small is "very small" in fact? For instance, should we take into account 
the contribution of 72 to jet cross-section in the region > 5 GeV, P 2 > 1 GeV 2 , where 
most of the data on virtual photons extracted from ep collisions at HERA come from? Despite 
the fact that the structure of longitudinal photons has not yet been investigated experimentally 
and is also poorly known theoretically, simple parton model estimates (see for details) of 
their effects suggest that in the mentioned kinematical region 7^ must be taken into account in 
the resolved photon contribution but may be safely neglected in the direct one. This difference 
comes from the fact that at small P 2 the contributions of to physical cross-sections behave 
as P 2 /s (i.e. vanish for fixed P 2 when s — > 00) in the direct channel, but as P 2 /^ 2 (with fi a 
fixed parameter) in the resolved part. In simple QED based considerations /1 is given by quark 
masses, while in realistic QCD we expect it to be given by some nonperturbative parameter 
of the order of 1 GeV. To assess the importance of 72, LO expressions for D e s(x, P 2 , M 2 ) 
evaluated with SaSID parametrizations for 7^, are compared in Fig. ^| with the formula (44) 
of for the quark distribution function of 72, treating m in the latter as a free parameter. 
The dotted curves in this figure correspond (from below) to m 2 = 1 GeV 2 , m 2 = 0.1 GeV 2 and 
to the asymptotic expression (a/2ir)l2e 2 x(l — x), obtained in the limit m — > 0. As expected, 
the importance of 72 depends sensitively on m. Moreover, its contributions relative to those of 
7^ peak at about x ~ 0.65, drop with increasing M 2 (for P 2 and m 2 fixed), and with increasing 
m 2 (for P 2 and M 2 fixed). Fig. ^ suggests that in the kinematical region accessible at HERA 
the contributions of 72 should be taken into account in phenomenological analyses of data 
on virtual photon interactions. Better theoretical knowledge of the structure of 72 is clearly 
needed. 

In all considerations of this subsection we simply added the contributions of 7^ and 72, 
despite the fact that the respective fluxes f^(y,P 2 ), f2(y,P 2 ) differ. Our conclusions are 
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Figure 6: D e s(x, P 2 , M 2 ) calculated from SaSID parametrizations for 7^ (solid curves) com- 
pared to results for^*L displayed by dotted curves and corresponding from above to m 2 = 0,0.1 
and 1 GeV 2 . 

therefore directly relevant for small values of y only but it is trivial to modify the above 
considerations by taking the respective fluxes properly into account. 



4 PDF of the virtual photon in NLO QCD calculations 



In [|J data on dijet production in the region of photon virtualities 1 < P 2 < 50 GeV 2 , and for 
jet transverse energies E^ > 5 GeV have been analyzed within the framework of effective PDF 
defined in (]T^) . This analysis shows that in the kinematical range 1 GeV 2 <P 2 <E 2 the data 
agree reasonably with the expectations based on SaS parametrizations of PDF of the virtual 
photon. The same data may, however, be also analyzed using the NLO parton level Monte- 
Carlo programs that do not introduce the concept of PDF of virtual photons, like DISENT, 
MEPJET or DISASTER++. Nevertheless, so long as P 2 <C M 2 ~ E\, the pointlike parts of 
PDF incorporate numerically important effects of a part of higher order corrections, namely 
those coming from collinear emission of partons in Fig. [I]. To illustrate this point we shall now 
discuss dijet cross-sections calculated by means of JETVIP currently the only NLO parton 
level Monte-Carlo program that includes both the direct and resolved photon contributions. 

All the above mentioned parton level NLO MC programs contain the same full set of partonic 
cross-sections Q for the direct photon contribution up the order act 2 . In addition JETVIP 

3 Note, however, that even at y = 0.4 the ratio f^,(y,P 2 )/f2(y,P 2 ) = 1-16 is still quite close to 1. 
4 In this subsection the various terms considered are characterized by the powers of a and a s appearing 
in hard scattering cross-sections. Writing O(a^a^) will thus mean terms proportional to aPa^ not terms up 
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includes also the resolved photon contribution. Once the concept of virtual photon structure is 
introduced, part of the direct photon contribution, namely the splitting term (H), is subtracted 
from it and included in PDF appearing in the resolved photon contribution. To avoid confusion 
we shall use the term "direct unsubtracted" (DIR uns ) to denote complete NLO direct photon 
contribution and reserve the term "direct" (DIR) for the results after this subtraction. In this 
terminology the complete NLO calculations is given by the sum of direct and resolved parts 
and denoted DIR+RES. 

For complete 0(a 2 ) calculations only the LO resolved photon contribution must be added to 
the 0{aa 2 s ) direct one. However, for reasons outlined in the next subsection, JETVIP includes 
also NLO resolved contributions. This might seem inconsistent as the corresponding 0{aa z s ) 
direct photon terms are not yet available and thus not included. Nevertheless, this procedure 
makes sense precisely because of a clear physical meaning of PDF of the virtual photon! 

4.1 Factorization mechanism in 7p interactions 

The main argument for adding 0(af) resolved photon terms to 0(aa 2 ) direct and 0(a 2 ) re- 
solved photon contributions is based on specific way factorization mechanism works for processes 
involving initial photons. This point is crucial but subtler and we therefore merely summarize 
the conclusions and refer the reader to the website version of this contribution, or to 0, for 
details. In the absence of 0(aa?) direct calculations, we have two options: 

• To stay within the framework of complete 0(a 2 ) calculations, including the LO resolved 
and NLO direct contributions, but with no mechanism for the cancellation of the depen- 
dence of PDF of the virtual photon on the factorization scale M. 

• To add to the previous framework the 0{a i s ) resolved photon contributions, which provide 
necessary cancellation mechanism with respect to the part of factorization scale depen- 
dence of photonic PDF "generated" by the homogeneous part of the evolution equations 
(jl]-|3|). The drawback of this procedure is the fact that the NLO resolved photon terms 
do not represent a complete set of 0(al) contributions. 

In our view the second strategy, adopted in JETVIP, is more appropriate. In fact one can look 
at 0(a^) resolved photon terms as results of approximate evaluation of the so far uncalculated 
0(aa z s ) direct photon diagrams in the collinear kinematics. There are of course 0{aa z s ) direct 
photon contributions that cannot be obtained in this way, but we are convinced that it makes 
sense to build phenomenology on this framework. 

Note that for the 0(a 2 ) resolved terms the so far unknown 0(aal) direct photon contri- 
butions provide the first chance to generate pointlike gluons inside the photon. To get the 
gluon in O(a^) resolved photon contributions would require evaluating diagrams of even higher 
order 0(aa^)\ In summary, although the pointlike parts of quark and gluon distribution func- 
tions of the virtual photon are in principle included in higher order perturbative corrections 
and can therefore be considered as expressions of "interactions" rather than "structure" , their 
uniqueness and phenomenological usefulness definitely warrant their introduction as well as 
their names. 

to this order! For approximations taking into account the first two or three powers of a s , in either direct or 
resolved channel, the denomination NLO and NNLO are used. 



9 



4.2 Dijet production at HERA 



To make our conclusions potentially relevant for analyses of HERA data on jet produnction 
we have chosen the following kinematical region (jets with highest and second highest are 
labelled "1" and "2") 

Jet transverse energies: asymmetric cuts — 7 GeV, E^ > 5 GeV; 

Photon virtuality: four windows in P 2 

1.4 < P 2 < 2.4 GeV 2 ; 2.4 < P 2 < 4.4 GeV 2 ; 4.4 < P 2 < 10 GeV 2 ; 10 < P 2 < 25 GeV 2 



Jet pseudorapidities in 7* CMS: -2.5 < 77® < 0, 2 = 1,2 



The cuts were chosen in such a way that throughout the region P 2 <C Ej., thereby ensuring that 
the virtual photon lives long enough for its "structure" to develop before the hard scattering 
takes place. As far as jet transverse energies are concerned, we have chosen the asymmetric cut 
scenario: E^r > E^ + A, EjP > E^, which avoids the problems || coming from the region 
where £jf ~ Ejt\ The asymmetric cut option is appropriate if one plots separately the E$p 
and EjP distributions, or, as we do, their sum, called distribution of "trigger jets" 0. To 
determine the value of A optimally, we evaluated the integral cr(A) over the selected region in 
Efp — E^p plane as a function of A and on the basis thereof took A = 2 GeV for all P 2 . 

In our analysis jets are defined by means of the cone algorithm. At NLO parton level all jet 
algorithms are essentially equivalent to the cone one, supplemented with the parameter -R S ep, 
introduced in order to bridge the gap between the application of the cone algorithm to NLO 
parton level calculations and to hadronic systems (from data or MC), where one encounters 
ambiguities concerning the seed selection and jet merging. In a general cone algorithm two 
objects (partons, hadrons or calorimetric cells) belong to a jet if they are within the distance 
R from the jet center. Their relative distance satisfies, however, a weaker condition 

AR V = ^A^+fA^ < J*IX) R - (U) 

The parameter R sep governs the maximal distance between two partons within a single jet, i.e. 
two partons form a jet only if their relative distance ARij satisfies the condition 



ARij < min 



Et ~\~ Et 

IX , -Tlsep 



max(E n ,E T 



(12) 



The question which value of -R sep to choose for the comparison of NLO parton level calculations 
with the results of the cone algorithm applied at the hadron level is nontrivial and we shall 
therefore present JETVIP results for both extreme choices R sep = R and R sep = 2R. To define 
momenta of jets JETVIP uses the standard ^-weighting recombination procedure, which leads 
to massless jets. To assess the reliability of our conclusions we have furthermore investigated 
the following uncertainties: 
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Choice of PDF: We have taken CTEQ4M and SAS1D sets of PDF of the proton and photon 
respectively as our canonical choice. Both of these sets treat quarks, including c and b 
ones, as massless above their respective mass thresholds, as required by JETVIP, which 
uses LO and NLO matrix elements of massless partons. Taking into account our cuts on 
jet transverse energies, we set Nf = 4 in any calculations discussed below. PDF of the 
proton are fairly well determined from global analyses of CTEQ and MRS groups and 
we have therefore estimated the residual uncertaintly related to the choice of PDF of the 
proton by comparing the CTEQ4M results to those obtained with MRS(2R) set. The 
differences are very small, between 1.5% at 77 = and 3% at 77 = —2.5. 

Factorization scale dependence: In principle proton and (in resolved channel) photon are 
associated with different factorization scales M p and M 7 , but we followed the standard 
practice and set M p = M 7 — M — k(E^ + E T 2 ')/2. The factorization scale dependence 
was estimated by performing the calculations for k = 1/2, 1, 2. 

Re normalization scale dependence: The dependence of perturbative calculations on the 
renormalization scale /1 is in principle a separate ambiguity, unrelated to that of factor- 
ization scale M, but we followed the common practice and set fi = M. 

Jet algorithm ambiguities: To see how much our conclusions depend on R sep we performed 
our calculations for two extremes values: R sep = R, 2R. 

Hadronization corrections: For any comparison of parton level calculations with exper- 
imental data understanding these corrections is crucial, but the problem is particularly 
pressing for jets with moderate Et, like those used in all analyses of virtual photon struc- 
ture. Hadronization corrections are not simple to define, but adopting the definition used 
by experimentalists , we have found that they depended sensitively and in a correlated 
manner on transverse energies and pseudorapidities of jets. For E T = 5 GeV they start 
to rise steeply below rj = —2.5 and we therefore required both jets to lie in the region 
—2.5 < 77 W < 0, where hadronization corrections are flat in 77 and do not exceed 20%. 



For virtual photons JETVIP can be run in two different modes: 



DIR uns : only the NLO unsubtracted direct photon calculations are performed without intro- 
ducing the concept of virtual photon structure. 

DIR+RES: employs the concept of PDF of the virtual photon and gives jet cross-sections 
as sums DIR+RES of subtracted direct and resolved photon contributions. 

We first discuss results for the first window 1.4 < P 2 < 2.4 GeV 2 . In Fig. [7|dcr/d77 and dcr/d-Er 
distributions of trigger jets obtained within the DIR+RES approach are compared to those of 
DIR uns one. All curves correspond to R sep = 2. The difference between the results of these 
two approaches is significant in the whole range of 77, but particularly large close to the upper 
edge 77 = 0, where the DIR+RES results exceed the DIR uns ones by a factor of more than 3! In 
dcx/dEV distributions this difference comes predominantly from the region of Et close to the 
lower cut-off E T + A = 7 GeV. Fig. [j] also shows that the scale dependence is nonnegligible 
in both approaches, but does not invalidate the main conclusion drawn from this comparison. 
The dependence of the above results on R sep (not shown) is almost imperceptible for DIR uns 
calculations and below 10% for the DIR+RES ones. To track down the origins of the observed 
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Figure 7: Scale dependence of the distributions da/dr] and der/cLEr o,t the NLO. All curves 
correspond to R scp = 2R. The dotted and dashed-dotted curves have the same meaning for 
DIR+RES as well as DIR uns calculations. 




Figure 8: Comparison of DIR+RES, DIR mis and DIR results for da/dr] (a). In b) individual 
contributions to da^l Q /dr] are plotted. 
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large differences between DIR+RES and DIR uns results, we did two simple exercises. In Fig. 
|]a the DIR+RES and DIR uns results are compared to the subtracted direct (DIR) ones. The 
difference between the DIR+RES and DIR curves, defining the resolved photon contribution 
da res /df], is then split into the contributions of: 

• the VDM part of photonic PDF convoluted with complete NLO (i.e. up to the order 
O(a^)) parton level cross-sections (denoted NLO VDM), 

• the pointlike quark and gluon distribution functions convoluted with 0(a 2 ) and O(a^) 
parton level cross-sections 

and plotted in Fig. [8|b. Fractional contributions of LO (i.e. 0(a 2 )) and NLO (i.e. 0(a 2 ) + 
0(o? s )) terms to <7 res are plotted in Fig. ^a as functions of r\. Several conclusions can be drawn 
from Figs. 

• The contribution of the VDM part of photonic PDF is very small and perceptible only 
close to r] = 0. Integrally it amounts to about 3%. Using SaS2D parametrizations would 
roughly double this number. 

• The inclusion of O(a^) resolved photon contributions is numerically important in the 
whole range —2.5 < rj < 0. Interestingly, throughout this interval the 0(a 3 ) results, 
particularly those of pointlike quarks, are close to the 0(a 2 ) ones. 

• At both 0(a 2 ) and O(ct^) orders pointlike quarks dominate da res /drj at large negative rj, 
whereas as r] — > the fraction of d<T res /d?7 coming from pointlike gluons increases towards 
about 40% at r] = 0. 



We reiterate that pointlike gluons carry nontrivial information already in convolutions with 
0(a 2 ) partonic cross-sections because in unsubtracted direct calculations such contributions 
appear first at the so far uncalculated order aa 3 . The results of pointlike quarks convoluted with 
0(«g) partonic cross-sections would be included in unsubtracted direct calculations starting at 
the order aa^, whereas for pointlike gluons this would require evaluation of unsubtracted direct 
terms of even higher order aaf! In JETVIP the nontrivial aspects of taking into account the 
resolved photon contributions can be characterized [] by the "nontriviality fractions" R3, R4 

_ q PL ® < s (0(c^)) + G PL ® a^(Q(a 2 )) _ G PL ag s (Q(^ 3 )) 
rt3 = , -K4 = , ( lo) 

^rcs (J res 



plotted as functions of 7] and P 2 in Fig. |K| Note that at rj = almost 70% of <7 rcs comes from 
these origins. This fraction rises even further in the region rj > 0. 

So far we have discussed the situation in the first window 1.4 < P 2 < 2.4 GeV 2 . As P 2 
increases the patterns of scale and R sep dependences change very little. On the other hand, the 
fractions plotted in Fig. |] and [10] vary noticably: 

• The unsubtracted direct photon contributions (DIR uns ) represent increasing fractions of 
the full NLO JETVIP results. 



'Disregarding the VDM part of resolved contribution which is tiny in our region of photon virtualities. 
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Tl Tl Tl 

Figure 9: Fractional contributions to a res . Upper and lower dotted (dashed-dotted) curves 
correspond to pointlike quarks (gluons) convoluted with NLO and LO partonic cross -sections. 
The lower dashed curve in a) gives the NLO VDM contribution. The solid curves show the 
ratia of NLO DIR uns calculations to the NLO DIR+RES ones. 




-2.5 -2 -1.5 -1 -0.5 -2.5 -2 -1.5 -1 -0.5 

Figure 10: Nontriviality fractions R 3 and R 4 as functions of rj and P 2 . 

• The relative contribution of pointlike gluons with respect to pointlike quarks decreases. 

• The nontriviality factor R4 (which comes entirely from pointlike gluons) decreases, whereas 
i?3, which is dominated by pointlike quarks and flat in 77, is almost independent of P 2 . 

All these features of JETVIP calculations reflect the fundamental fact that as P 2 rises towards 
the factorizations scale M 2 « the higher order effects incorporated in pointlike parts of 
photonic PDF vanish and consequently the unsubtracted direct results approach the DIR+RES 
ones. The crucial point is that for pointlike quarks and gluons this approach is governed by 
the ratio of P 2 /M 2 . The nontrivial effects included in PDF of the virtual photon thus persist 
for arbitrarily large P 2 , provided we stay in the region where P 2 <C M 2 . 
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5 Summary and conclusions 



We have analyzed the content of parton distribution functions of the virtual photon within 
the framework of SaS approach to their parametrization. In this approach, quark as well as 
the gluon distribution functions can be separated into the nonperturbative VDM and pointlike 
parts, the latter being in principle calculable by perturbative means. The inherent ambiguity of 
this separation, numerically large for the real photon, becomes phenomenologically irrelevant 
for virtual photons with P 2 > 3 GeV 2 . In this region quark and gluon distribution functions 
of the virtual photon are dominated by their (reasonably unique) pointlike parts, which have 
clear physical origins. We have analyzed the nontrivial aspects of these pointlike distribution 
functions and, in particular, pointed out the role of pointlike gluons in leading order calculations 
of jet cross-section at HERA. 

At the NLO we have found a significant difference between JETVIP results in approaches 
with and without the concept of virtual photon structure. While for the real photon analogous 
difference is in part ascribed to the VDM components of photonic PDF, for moderately virtual 
photons it comes mostly from pointlike quarks and pointlike gluons. Although their contri- 
butions are in principle contained in higher order direct photon calculations, to get them in 
practice would require calculating at least 0(aa^) and 0(aaf) unsubtracted direct contribu- 
tions. In the absence of such calculations the concept of PDF of the virtual photon is therefore 
phenomenologically indispensable. 

We have also shown that despite the expected dominance of the transversely polarized 
photons, the longitudinal polarization of virtual photons may play a numerically important 
role in the region of moderate virtualities accessible at HERA and should therefore be taken 
into account in analyses of relevant experimental data. 

Acknowledgments: We are grateful to G. Kramer and B. Potter for discussions on the 
interactions of virtual photons and to B. Potter for help with running JETVIP. 
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